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Abstract 

Proposed in this paper is a numerical procedure to generate periodic traveling wave 
solutions of some nonlinear dispersive wave equations. The method is based on a 
suitable modification of a fixed point algorithm of Petviahvili type and solves several 
drawbacks of some previous algorithms proposed in the literature. The method is 
illustrated with the numerical generation of periodic traveling waves of fractional 
KdV type equations and some extended Boussinesq systems. 
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1 Introduction 


Considered here is a technique to generate numerically periodic traveling wave 
solutions of nonlinear wave equations of the form 


ut - Lux + f{u)x = 0 , X eR,t>0. ( 1 ) 

In (II]) M = u{x,t), / : M —)■ M is a smooth function (typically a polynomial) 
with /(O) = f'{0) = 0 and L is a dispersion operator with Fourier symbol 
a : M —)■ M, 


LviO = aiom, e e M, (2) 

being v{^) the Fourier transform of n at G M. Equation ([T]) appears as a 
hrst one-dimensional, unidirectional approach for nonlinear dispersive wave 
propagation, see e. g. 011 and particular cases are remarkable: Korteweg 
de Vries-type equation (KdV, L = —dxx^fiu) = > 3), Benjamin-Ono- 

type equation (BO, L = —l-idx with Ti the Hilbert transform, f{u) = > 

3) or Benjamin-type equation (BO, L = —jSdxx — I'Hdx with > 0, f{u) = 
uP~^,p > 3). 

The paper is concerned with the numerical generation of periodic traveling 
wave solutions of the models ([T]). They are solutions of the form 


u{x,t) = (j){x - Cst) = (3) 

n=—oo 

whic are periodic of some period 21,1 > 0 and travel with speed Cs > 0. 
Results on existence and stability of such solutions for particular cases of ([T]) 
can be seen in the literature (see e. g. [6] and references therein). Recently, 
[12], Chen and Bona extablished existence of periodic traveling wave solutions 
of (JI]) under some hypotheses on / and L (see below) as well as stability under 
perturbations of the same period. This study generalizes previous works on 
this matter, [12]. 

In [2] a numerical scheme for the numerical generation of periodic travel¬ 
ing waves of (II]) was proposed. This was based on a Fourier pseudospectral 
discretization in space of ([T]) along with time integration with conservation 
properties of invariant quantities of the periodic value problem associated to 
dl]). The periodic traveling wave prohle was approximated by using the iter¬ 
ative Petviashvili method, mm- The resulting scheme ensured an accurate 
computation of two characteristic elements of the waves, the initial prohle 
and the speed. The present work attempts to complete [2] in two main points. 
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not totally resolved there in the authors’ opinion. The first one is that the 
generation of the prohle of the periodic traveling wave solution may involve 
some nonzero constants which make the Petviashvili method, at least initially, 
unable to be used. The second one, somehow related, is that this method is 
not well suited to nonlinearities / in ([1]) with more than one homogeneous 
terms of different degrees (like a polynomial for example). 

In order to overcome these drawbacks, the present work introduces an im¬ 
proved algorithm from the one proposed in [2]. The strategy is based on a 
change of variable suggested in [13] (see also m) to treat the constants (to 
solve the hrst problem above described) and the use of extended versions of 
the Petviashvili method introduced in |3] that allows the presence of more 
general nonlinearities in the model (to solve the second problem). The struc¬ 
ture is as follows. In Section |2] some hypotheses on ([ID (considered in [13]) are 
assumed and the proposed algorithm is decribed. The efficiency of the method 
is checked in Section [3]by considering the periodic traveling wave computation 
for two models. The hrst one is the generalized fractional KdV equation (of 
the form ([ID with Q!(,^) = < /r < 2,/(m) = — 1)), [TIITS] . The 

second model is the e-Boussinesq system, [13]. This is not of the form ([ID but 
its presence here is justihed as an attempt to extend the application of the 
method is based to more general situations, in particular to some nonlinear 
dispersive systems. In both cases and to our knowledge, many of the prohles 
are generated for the hrst time; in particular, for the case of the e-Boussinesq 
system, we do not even know if some theoretical results of existence of periodic 
traveling wave solutions for the model are available (although they may be ex¬ 
pected from the existence of solitary wave solutions as their large wavelength 
limit, i)- Section 0] contains some concluding remarks. 


2 Preliminaries and numerical method 


We make the following assumptions about ([ID: 

(HI) The nonlinear term / is a polynomial of the form 

f{z) = 373^2 -!-•••+ pipzP-^, 

where 7j > 0, j = 3, ... p — 1, 7p > 0, p > 3. 

(H2) a is real, even, continuous with q;(0) = 0. If po = min{j = 3,... ,p/'^j > 0} 

_2 

then there exists s >-such that 


lim ^ = 0. 

|^| 2 . 
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(H3) There exists s > — such that 



These are some of the hypotheses considered in |T3] and will be adopted here 
with the aim of ensuring the existence of periodic traveling wave solutions in 
the following sense (see |T3] for details): For a sufficiently large period param¬ 
eter I > 0 the prohles 0 are obtained as minimizers of constrained variational 
problems involving momentum and energy functionals for ([T]). Then for each 
minimizer 0 one can hnd a positive speed Cg such that ([3]) is an inhnitely 
smooth traveling wave solution of ([T]) which is also stable under perturbations 
of the same period. In [13] the authors also study the large wavelength limit 
/ —)■ oo of the periodic traveling waves determining under some additional 
hypotheses the convergence to solitary wave solutions. 

A comment in [13] is the starting point to introduce our numerical strategy 
to generate numerically periodic traveling wave prohles. When n in (I3|) is a 
solution of ([T]) then the prohle 0 must satisfy 


c^0-L0 + /(0) = A, 


(4) 


where A is a (nonzero in general) constant. Then our method is based on the 
following steps: 

(SI) Finding (approximate or exact) real constants solutions 0 = C* of (jl]) that 


IS 


+ f(C) = A. 


(5) 


(Note that LC = 0 since q:(0) = 0 is assumed in (H2).) This means that C 
must be a root of the polynomial 


P(z) = -A - CgZ + Sjsz'^ H-h pjpZ^ ^ 


(S2) The change of variables (p = (j) — C, when 0 is a solution of (jl]) along with 
([5]) lead to the following equation for ip-. 



L{C) 


(6) 


N(f) 


where 
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p 




1=3 


P 


f\C) = j;: !(! - 1) ■ ■. (i - i)'nC'-=-\ j > 2. 


l=j+l 


(S3) Equation ([H]) is iteratively solved with a extended version of the Petviashvili 
method, [1]. Writing 



and from ^ 0 and for z/ = 0,1,... the following iteration is considered: 


p-i 


£( C ) 4,[‘'+‘1 = = 0 , 1 ,..., 


(7) 




The factors Sj in (0) are obtained by using the same quotient s{u) but 
changing the exponent aj. Formulas (I7j), (IH]) correspond to an extension of 
the classical Petviashvili scheme for nonlinearities with several homogeneous 
terms but different degree of homogeneity. This belongs to a more general 
family, described in |1]. The use of O. B instead of other alternatives 
presented in the literature (of Newton-type, for example) is justihed by the 
attempt to simplify the implementation of the iterative techniques for the 
prohles with an easily implementable hxed point type algorithm, [19] . Note 
that in the case of convergence, —?• 1. 

The operator ([2D taking part of L{C) in ([6D suggests to consider the imple¬ 
mentation of (ED, ([HD in Fourier space for the Fourier coefficients G 'L 

of the iterations <^9^, z/ = 0,1,... Then for hxed Cg and I > 0 the problem ([ID 
on {—1,1) with periodic boundary conditions is considered and (ED, ([HD fake 
the form 



1 


p-i 


(9) 


n G Z, z/ = 0,1,..., 



(10) 


j = 2,...,p-l, 
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Finally, we note that the existence results in [13] are mentioned here in connec¬ 
tion with the hypotheses (H1)-(H3). However, the strategy for the numerical 
generation uses a different argument, in the sense that it previously hxes a 
speed Cg whose associated prohle is attempted to be computed. In this sense, 
the point of view of the resolution is closer to other existence results presented 
in the literature, |7]. 


3 Numerical examples 


This section is devoted to illustrate the behaviour of the algorithm to generate 
numerically periodic traveling waves. This will be shown by considering two 
models: a generalized version of the fractional KdV equation and the extended 
Boussinesq system. 

It may be worth making hrst some brief comments on the implementation of 
O, dHI). As mentioned above, due to the nonlocal operator (]2D, a natural way 
to approximate ([6D consists of a discretization on an interval with a Fourier 
pseudospectral method as in [2], see [10]. This allows to implement the iter¬ 
ation by using the discrete Fourier version of ([9]), flTOD and FFT techniques. 
Secondly, in some of the experiments, the resulting scheme has been com¬ 
plemented with vector extrapolation methods to accelerate the convergence. 
The literature on these methods is abundant (see e. g. mm and references 
therein) and their use in traveling wave simulations is analyzed in [5]. 


3.1 Generalized fractional KdV equation 

The hrst example treated here concerns the generation of periodic traveling 
waves of fractional KdV models of the form 


ut - (/(m))x = 0, 


( 11 ) 


which corresponds to ([T]) with L = A^, being 


Ar^(0 = 


u 


p-i 


0 < /i < 2 and f(u) = - ,p > 3. Existence and stability 

p — 1 

traveling wave solutions are analyzed in [TB] by using variational 


of periodic 
arguments. 
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Specifically, existence is obtained for 0 < p < Pmax where 


p < 1 

/r > 1 

As far as stability is concerned, for /i > 1/2 and 

P =-2+(^-1)2.-' 

then the periodic traveling wave is 

• Spectrally stable if p > 1 and 1 < p < p*(p). 

• Spectrally unstable ifl/2</i<lorp>l when p > p*(p). 

Furthermore, if 

P*max = maxp*(p), 

/i>l 

then for p G (l,Pmax) there are /i-(p) < p+(p) such that small amplitude peri¬ 
odic traveling waves are spectrally stable for p G (p_(p), p+(p)) and spectrally 
unstable for a ^ (p_(p),p+(p)). 

In this case step [SI] of the algorithm requires the computation of constant 
solutions u = C satisfying 


Pmax 



-l-CXD if 


CP-i 

- CsC -A = 0, A G M. (12) 

p — 1 

Existence of real solutions of flT^ requires some relations involving the pa¬ 
rameters p^Cs^A (but not p). From the constant C of flT^ the equation ([6]) 
takes the form 


(c, + L-C'P-2)p 

V 

L(C) 



= 0 , 


N{^) 


The behaviour of the algorithm is illustrated here by hxing A = = 1 and 

two values for a and p\ a = 0.8, p = 3 and a = 1.5, p = 4. In both cases, the 
convergence was accelerated by using the minimal polynomial extrapolation 
method (MPE), [TT1IT8] . The approximate prohle along with the corresponding 
phase portrait, are shown in Figures [T](a),(b) and[2][a),(b), respectively. Figures 
[T][c),(d) and [2](c),(d) check the accuracy by displaying two types of errors as 
functions of the number of iterations: 
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(a) 


(b) 



Fig. 1. Periodic traveling wave generation of CID with A = Cs = 1 and 
p = 3, a = 0.8. (a) Approximate profiles; (b) Phase portrait; c) Discrepancy (fT4j) 
vs number of iterations; (d) Residual error (|13p vs number of iterations. (Semi-log¬ 
arithm scale in both cases.) 


• The residual error (in Euclidean norm) 

RES, = IlLypH - iV((^M)||^ z/ = 0,1,... (13) 

• The discrepancy between the sequence of stabilizing factors s, = and 

(in case of convergence) its limit one 

SFE, = \su-l\, z/ = 0,1,... (14) 


3.2 Extended Boussinesq system 


The second example attempts the application of the method to systems of the 
form ([T]). Specihcally we consider the so-called e-Boussinesq system 








(a) 


(b) 



Fig. 2. Periodic traveling wave generation of m with A = Cs = 1 and 
p = 4, a = 1.2. (a) Approximate profiles; (b) Phase portrait; c) Discrepancy (fT411 
vs number of iterations; (d) Residual error (I13p vs number of iterations. (Semi-log¬ 
arithm scale in both cases.) 

r]t = -diW^ - d2Wxxx - d4{Wr])x + d^{Wrf)x, (15) 

Wt = - dsWxxt - + d,{W^v)., (16) 


Equations ffTHll . ffTHll appear as a model for the bidirectional propagation of 
an interfacial wave ri{x,t) between two fluid layers with W{x,t) standing for 
the horizontal velocity of the flow. The constants are defined as 


di 

ds 


H 


r + H' 




sdi ^ ^ _ r{l + H)^ 

“ (r + hf)2’ ~ (r + 


where r and H are, respectively, the (dimensionless) density and depth ratios 
while s is related to the physical depth below and above the unperturbed in¬ 
terface and is taken within the range —{1 + rH) < s < —^{1 + rH). The model 
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is described in detail in [16], where the propagation and collision of solitary 
waves and fronts are also studied numerically. To our knowledge existence of 
periodic traveling wave solutions has not been analyzed (theoretically or nu¬ 
merically). Suggested here is some numerical evidence of such waves by using 
the method described in Section |2l 

Note hrst that periodic traveling waves 7] = r]{X), W = W (X), X = x — Cgt of 
firsj) . ffT6|l of speed Cg must satisfy the system for the prohles 


Cg (di -|- d2dxx)\ ^ 


Cg{\^d^dxx) ) \W) \-W\-^ + d5v) 


-Wrj^-di d^rj) 


+ 




(17) 


for some integration constants 741,^42. The search for constant solutions r] = 
Cl, VT = C 2 of (ITTII leads to the algebraic system 


^ Cg -di'^ ^ 


\-T, 


Cl 


^d4CiC2 - d^ClC2\ ( Ai 

^fCl-dsWlj V7l2 


in such a way that if 7 ^ 1 then 


C 2 = 


CgCi — Ai 


di + d4Ci - dsC'r 
while Cl must be a root of the polynomial 


P{z) = ——z{di + d^z — d^z^Y -T Cg{cgZ — Ai){di + d^z — d^zY 


di 


d/i 


2\2 


~\~{,d^z — ^)(cs^ ~ ^1) ~ ^2(^1 A diZ — d^z ) 


Finally, the system for the differences rj = rj — Ci,W = IF — C 2 is of the form 


\-d^ + 2d5Ci)C2 CY-dA + d,CY 

^ d,Cl 2 C 2 (-f+ d5C'i) 

d5C2fY + Wrj{-di + d^rj) + d^WrY 

^lF2(-f + d^fj) + 2d^C2Wfj + dYA^^rj 



Cg —{di + d2dxx) 
+Cs(l + d^dxx) 


\wj 

(18) 


Then, as mentioned in Section [J] system (ITSU is iteratively solved for rj, W and 
approximate periodic traveling wave prohles for the values Ai, A 2 are obtained 
from the formulas rj = fj + Ci,W = hF-|-C 2 . 
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In [16] a maximum wave velocity 


- 1 + + 

is derived and solitary waves are characterized by wave velocities larger than 
one. We illustrate the performance of the method by generating periodic trav¬ 
eling waves with a speed close to Vmax (specihcally Cg = Vmax — 10““^ in the 
experiments below) and two different values for r, H and the constants ^42. 
(In both experiments s = —(1 -I- rH) has been hxed.) The hrst case, Figure 



X T1 


(a) (b) 




Fig. 3. Periodic traveling wave generation of (1T71) with Cg = Vmax — 10“'^. 
Ai = —1,^2 = —2,r = 0.8, H = 0.95, s = —1 — rH. Vmax = 1.000454. (a) Approxi¬ 
mate profiles; (b) rj phase portrait; (c) Residual error (fT3ll vs number of iterations; 
(d) Discrepancy (I13|) vs number of iterations (semi-logarithmic scale). 

[3l corresponds to r = 0.8, FT = 0.95, Ai = —1,^2 = —2. Figure [3]( a) shows 
the approximate rj, W prohles (also obtained by using additional acceleration 
with the MPE method) and the periodic character is conhrmed by the phase 
portrait of r] in Figure |3](b). The accuracy of the computation is again evalu¬ 
ated in Figures |3](c) and [3](d) with the behaviour of the errors flTOll and flTTl) 
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with respect to the number of iterations. 


The same information is provided in Figure 0] for the case r = 0.8, H = 
1.8, y4i = 1,^2 = 1. Note here that the periodic form of the approximate 
Tj prohle is different and closer to the ‘table-top’ form of the solitary waves 
shown in [16]. This suggests the connection between these periodic prohles 
and some associated solitary wave in the limit of large wavelength, [9]. 
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(a) 


(b) 



Fig. 4. Periodic traveling wave generation of (I17p . Approximate profiles and phase 
portraits with Cg = Vmax — 10“^. Ai = 1, A 2 = l,r = 0.8, H = 1.8, s = —1 — rH. 
Vmax = 1.065919. (a) Approximate profiles; (b) rj phase portrait; (c) Residual error 
(jlOp vs number of iterations; (d) Discrepancy (I18p vs number of iterations (semi-log¬ 
arithmic scale). 


4 Concluding remarks 


In this paper an iterative technique to generate periodic traveling wave solu¬ 
tions of some nonlinear dispersive wave equations is proposed. The method 
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is based on a suitable change of variable to avoid the presence of some con¬ 
stants of integration in the problem and the application of hxed point iteration 
of Petviashvili type, |1]. This strategy overcomes two difficulties not totally 
solved in some previous algorithms, |2], for this purpose, namely: 

• The presence of constants that prevents the use of Petviashvili type meth¬ 
ods, which have been shown to be an alternative to compute solitary wave 
solutions for the models under study (see [3] and references therein). 

• The presence of nonlinearities of polynomial type with homogeneous terms 
of different degree, for which the construction of the Petviahsvili type meth¬ 
ods is not suitable. 

The scheme has been illustrated by computing periodic traveling waves in two 
problems of additional interest: the generalized fractional KdV equation and 
the extended Boussinesq system. In the hrst case, existence and stability of 
periodic traveling waves have been studied theoretically but, to the authors’ 
knowledge, no approximate prohles have been computed. The second exam¬ 
ple attempts to extend the application of the technique to some nonlinear 
dispersive systems and suggests the formation of periodic traveling waves in 
the model, for which (again to the authors’ knowledge) there are no existence 
results. 
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